A mean-field theory of Anderson localization 
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Anderson model of noninteracting disordered electrons is studied in high spatial dimensions. We 
find that off-diagonal one- and two-particle propagators behave as gaussian random variables w.r.t. 
momentum summations. With this simplification and with the electron-hole symmetry we reduce 
the parquet equations for two-particle irreducible vertices to a single algebraic equation for a local 
vertex. We find a disorder-driven bifurcation point in this equation signalling vanishing of diffusion 
and onset of Anderson localization. There is no bifurcation in d = 1, 2 where all states are localized. 
A natural order parameter for Anderson localization pops up in the construction. 
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Understanding mobility of (free) particles in random 
media has been a challenging theoretical problem for 
many decades. It became clear from the early days of 
the study of systems with randomly distributed scatter- 
ers that the particle movement has a diffusive charac- 
ter described at long distances by a diffusion equation. 
A breakthrough in the conceptual perception of random 
systems was achieved in Ref. [ij . P. W. Anderson demon- 
strated there on a simple model that for a sufficiently 
strong disorder the particle remains confined in a finite 
volume and fails to diffuse to long distances. Since then, 
the disorder-induced absence of diffusion, called Ander- 
son localization, has attracted a lot of attention of con- 
densed matter theorists. In spite of years of intensive 
studies, the phenomenon of Anderson localization has 
not yet been fully understood. It is mostly due to many 
facets of this difficult problem. 

There are two complementary tools, numerical and 
analytical, to tackle Anderson localization and the metal- 
insulator transition. Neither of them is, however, able 
to answer all questions about the disorder-induced van- 
ishing of diffusion. While the former deals with finite 
lattices and many configurations of the random poten- 
tial the latter deals mostly with the thermodynamic 
limit and configurationally averaged quantities [2|. Al- 
though the existence of Anderson localization was proved 
rigorously to exist in any finite dimension if the disor- 
der is sufficiently strong both numerical and ana- 
lytical approaches work preferably in rather low dimen- 
sions {d = 2 -\- e). As a consequence, a standard mean- 
field (high-dimensional) theory of Anderson localization 
is missing, except for special solutions on a Bethe lattice 
0. Even though a self-consistent, mean- field-type the- 
ory of Anderson localization was formulated |6|, it was 
derived only within the weak-scattering, low-dimensional 
limit. A systematic mean- field theory should come out 
of the asymptotic limit to high spatial dimensions ■ 

The aim of this Letter is to employ the limit to high 
spatial dimensions for developing a mean-field theory for 
the Anderson localization transition. We use the parquet 



approach summing up systematically nonlocal vertex cor- 
rections to the mean-field, d = oo, solution f§]. We show 
how the parquet equations for the irreducible two-particle 
vertices can be simplified and solved in the asymptotic 
limit to high spatial dimensions. We use the asymptotic 
solution of the parquet equations to build up a quantita- 
tive theory of the disorder-driven metal-insulator transi- 
tion with a mean-field critical behavior, i. e., independent 
of the spatial dimension. 

The existence or nonexistence of diffusion can be de- 
termined from the electron-hole correlation function de- 
fined from the two-particle resolvent as $(21,22; q) = 
A^^^ g|^^i^, (zi, Z2; q), where zi and Z2 are complex 
energies. A specific element of this function with ener- 
gies zi — Ep + uj + lO"*" and Z2 = Ep — zO"*", denoted 
$d^(q, cj), is used to determine the diffusion constant 



(1) 



where np is the density of states at the Fermi energy 
Ep. Vanishing of the diffusion constant D indicates the 
absence of diffusion in the system. 

In systematic theories we do not approximate directly 
either the Green function G^^-* or the correlation function 
<I>^^, but rather the two-particle vertex F defined from 



G+(k)G_(k + q) [(5(k-k') 

+Fkk'(q)G+(k')G_(k' + q)] 



where G±(k) = G(k, z±) are averaged one-particle resol- 
vents. 

The simplest theory for strongly disordered systems 
is the coherent-potential approximation (CPA), being an 
exact solution in rf = 00. Only the diagonal part of the 
one-electron resolvent G{z) = ^(^i ^) rele- 

vant in the CPA, since it completely neglects coherence 
between spatially distinct scatterings. It is a consistent 
mean- field theory only for one-electron functions. Nonlo- 
cal parts of two-particle functions in high dimensions do 
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not vanish. The ofF-diagonal elements of the one-electron 
propagators cannot be neglected and have to be taken 
expliciltly into consideration 7]. Further on, the CPA, 
due to its local character, is degenerate and cannot dis- 
tinguish between scatterings of electrons and holes. To 
be able to resolve various types of two-particle scatter- 
ings in noninteracting systems we have to go beyond the 
mean-field, local approximation. 

A systematic (diagrammatic) expansion around the 
d = oo limit can be performed by using the off-diagonal 
one-electron CPA resolvent G(k, 2) = G(k, z) — G(z). 
Three nonequivalent Bethe-Salpeter equations with non- 
local propagators can be constructed for the full vertex 
r 0. Here we will use only the Bethe-Salpeter equations 
from the electron-hole and the electron-electron scatter- 
ing channels They can be represented as 

rkk'(q)-Afk.(q) + ^E^kk"(q) 

^ (2a) 

X G+(k")G_(k" + q)rk"k'(q) . 
Fkk' (q) - A^k' (q) + ^ E Akl" (q + k' - k" ) 

k" 

X G+(k")G_(Q-k")rk"k'(q + k-k") , 

(2b) 

respectively. We used bar in the irreducible vertices A'^^ 
and A'^^ to indicate that the Bethe-Salpeter equations 
are constructed with the off-diagonal resolvents only. 
Hence, in the infinite-dimensional case G±(k) = and 
^e?i,ee _ where 7 is the full local CPA vertex 1^. For 
simplicity of notation we denoted Q = q-|-k + k'. Notice 
that q is the momentum conserved for scatterings in the 
electron- hole channel, Eq. (|2a|l . while Q is conserved in 
the electron-electron channel, Eq. (|2b(l . 

To reach the strong-disorder limit with a diffusionless 
regime we have to determine the irreducible vertices A'^^ 
and A*^*^ self-consistently. The parquet construction pro- 
vides a suitable framework for this purpose, ft is based 
on the observation that reducible diagrams in one channel 
are irreducible in the other, topologically distinct chan- 
nels. If we approximate the vertex irreducible in all chan- 
nels by the local CPA vertex 7, take into account only the 
eh and ee channels, and realize that the full vertex is a 
sum of reducible and irreducible vertices in any channel, 
we end up with a fundamental parquet equation 

Fkk'(q) = Afk'(q) + Ark'(q)-7. (3) 

The minus sign at the CPA vertex compensates for the 
identical local part in both A^^ and A^'^. A couple 
of (nonlinear) parquet equations determining the irre- 
ducible vertices A^^ and A*^"^ as functions of 7 and G± 
are obtained by replacing the full vertex F in Eqs. Q by 
Eq. ®. 

Prior to solving the parquet equations for A'^'* and 
A'^'^ we utilize the electron-hole symmetry expressed as 



an identity for two-particle vertices Fkk'(q) — rkk'(^q~ 
k - k') and A^^^, (q) = A^'^, (-q - k - k') . The electron- 
hole transformation maps Eq. (|2a|l onto Eq. (|2b(l . The 
two-particle electron-hole symmetry is a consequence of 
the time-reversal invariance the one-electron resolvent 
G(k, z) — G(— k, z) used in the Bethe-Salpeter equations 
(0) . This symmetry actually reduces the number of par- 
quet equations to a single nonlinear integral equation for 
Akk' (q) ^ A- , (q) = Ali, (-q - k - k') . 

Generally, the parquet equations are unsolvable due 
to momentum convolutions in the Bethe-Salpeter equa- 
tions. The limit to high spatial dimensions leads to sup- 
pression of spatial fluctuations resulting in simplifications 
of momentum convolutions We take advantage of 

these simplifications. We start with the leading asymp- 
totic term in the off-diagonal propagator G(k, z) that on 
a d-dimensional hypercubic lattice with the hopping am- 
plitude t reads 

4- r 
G(k, z)^^Yl cos(fc^) / dep{e)G{z ~ I](z) - ef 

= txik){G{zf) , (4) 

where p is the density of states and S the local (CPA) self- 
energy. We replace the off-diagonal one-electron prop- 
agators in the parquet equations with this asymptotic 
representation. 

The simplest and most important convolution is a 
two-particle bubble x(q) = A^"^ EkG'-h(k)G_(k -t- q). 
Its asymptotic behavior can be found from the following 
formula 

1 If'* 
- a;(k)x(k + q) - -A(q) = ^ E (^a) 

k u=l 

where we denoted X{q) a two-particle (bosonic) disper- 
sion function. Other possible convolutions of the generic 
fermionic and bosonic dispersion functions are 

^^A(q' + q)a:(q'+k) = ^x(q-k) , (5b) 
q' 

l^A(q' + qi)A(q' + q2) = ^A(qi-q2) . (5c) 
q' 

We can see that the fermionic and bosonic dispersion 
functions form a closed algebra with respect to momen- 
tum summations. The elementary convolutions Q man- 
ifest the generation of the factor d^^ due to mixing of 
two-particle propagations from different scattering chan- 
nels. The dispersion functions then behave in the lead- 
ing asymptotic order of c?~^ as gaussian random variables 
when momentum summations are performed. 

To make the calculations in high spatial dimensions 
more mean- field- like, we replace the bare fermionic and 
bosonic dispersion functions with the respective off- 
diagonal one- and two-particle propagators. That is, we 



3 



use G instead of x and x instead of X. These quan- 
tities are directly proportional in the leading asymp- 
totic order. Without loosing the asymptotic accuracy 
we can extend relations |(SJ) to genuine mean-field expres- 
sions Eq' X(q' + q)G±(q' + k) = WG±{ci ~ k)/4d, 
Eq X(q + qi)x(q + q2) = Wx{q.i - q2)/4d, where 
^eusedW = t^{G\){Gl_). 

It is evident from Eq. and Eq. lO that the two- 
particle vertices can be represented as functions of the 
generic off-diagonal fermionic G(k) and bosonic x{^) 
functions. To find the leading high-dimensional asymp- 
totics of the solution of the parquet equation for Akk' (q) 
we keep only the leading terms for each specific mo- 
mentum dependence of the vertex. It is easy to demon- 
strate that the parquet equation then simplifies in the 
leading asymptotic limit d ^ cx3 to an algebraic equation 



Akk'(q) = A(q) =7 + Ao 



Aox(q) 
1 - Aox(q) 



(6a) 



where Aq = Af^^EqA(q). The high-dimensional ir- 
reducible two-particle vertex is completely determined 
from a single local parameter Aq and the two-particle 
bubble x(q)- Summing both sides of Eq. (|6a|l over mo- 
menta we obtain 



Ao = 7 -I- Aq 



1 



Aox(q) 



TV ^ 1 - Aox(q) 



(6b) 



Equations ^ were derived from the leading high- 
dimensional asymptotics, but can be used in any finite 
dimension. It is, however, mandatory that the proper 
d-dimensional momentum summation is used on a d- 
dimcnsional lattice. We cannot directly use the gaus- 
sian rules to evaluate the summation over momenta in d 
dimensions as we did during the derivation of these equa- 
tions. The integrand would be singular (nonintegrable) in 
the gaussian evaluation. To assess the asymptotic high- 
dimensional behavior of the two-particle irreducible ver- 
tex we have to realize that in deriving Eq. © each inde- 
pendent one-dimensional momentum integration over the 
components qi with i = 1, 2, . . . , d can contain maximally 
squares of the dispersion to remain within the leading 
asymptotics. Hence, on a d-dimensional lattice we can 
build maximally d pairs of the two-particle bubbles x- 
The integrand in Eq. (|6b|) therefore collapses to a poly- 
nomial of order d. Using the gaussian integration rules 
we explicitly obtain 



Aox(q) 



1 

N^l- Aox(q) 



E 



(2n)! 
(2d)"n! 



(a) 



where we denoted a = VF^Ag/8. The asymptotic limit 
of the sum fd-,oo{a) converges for a < 1/2, i. e., for 
Aq < A/W'^. The critical value of the randomness 



7^ = A/W'^ defines an ultimate upper bound beyond 
which perturbation theory around d = cxd in powers of 
d^^ docs not converge and becomes nonanalytic. In real- 
istic models, however, such an extreme value cannot be 
reached, apart from tiny regions around band edges and 
in satellite bands. 

The one-electron functions G, S and the local two- 
particle vertex 7 entering Eqs. © were assumed to be 
taken from the CPA. There we have 7 = A/(l - AG+G_) 
and A = (I]_|_ — !]_)/ (G+ — G_ ) . If we analogously define 
Ao = Ao/(l - AoG+G_) and x(q) = x(q) + G+G_ we 
can represent the asymptotic form of the full two-particle 
vertex as follows 



kk 



(q) = 7 



An 



Aox(q) 



Aox(k + k' + q) 



Ll-Aox(q) l-Aox(k + k' 



(7) 



The full nonlocal CPA vertex can be recovered from the 
above expression if we put Aq = 7 and neglect the last 
term on the r.h.s. of Eq. The term neglected in 

the CPA, however, restores the electron-hole symmetry 
in the asymptotic vertex in high dimensions. 

Up to now we have analyzed the two-particle asymp- 
totics with the one-electron propagators fixed by the 
CPA. To reproduce diffusion in this approach we have to 
match correctly the irreducible vertex calculated from the 
parquet equations lO and the one-electron self-energy. 
We hence have to go beyond the CPA even in the one- 
electron propagators. To do so consistently we use the 
Ward identity and determine the imaginary part of the 
self-energy from the two-particle irreducible vertex via 

^Y.{E + iO+) =Ao{E + iO+, E - iO+)^G{E + iO+) . (8) 

The real part of the self-energy is determined from the 
Kramers-Kronig relation Isll . Equation © completes the 
parquet equation (j6b|) . Both equations together with 
the Kramers-Kronig relation have to be solved simulta- 
neously to achieve full self-consistence between Aq and G 
calculated from S via the Dyson equation. 

The redefinition of the self-energy in Eq. is impor- 
tant, since only with it we recover the diffusion pole in 
the vertex functions. That is, we obtain Ao(£'-l-iO+, E — 
iO"'")x(0) = 1, whenever the parquet equation H6b|l al- 
lows for a positive solution. The existence of the dif- 
fusion pole is essential for the diffusion constant from 
Eq. to be non-zero (positive). We can immediately 
conclude from simple power counting in the momentum 
integral of Eq. Ij6bp that there is no positive solution for 
Ao(£'-|-jO+, E—iO^) in low dimensions, d — 1,2, since the 
diffusion pole would be nonintegrable. Consequently, no 
diffusion pole can exist and the diffusion constant from 
Eq. Q vanishes in d = 1, 2. 

In higher dimensions we can expand the r.h.s. of 
Eq. f\6h\i in powers of the local vertex Aq. We then obtain 
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FIG. 1: DifFusion constant D and the order parameter in 
the locahzed phase SiAo calculated from Eq. @. We used a 
semielliptic energy band with the bandwidth 2u;, the self- 
consistent Born approximation for the self-energy, and set 
Cd = O.IW^ 
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FIG. 2: Phase diagram for the same setting as in Fig.Q The 
hatched area denotes localized states. 

a mean-field-like cubic equation 

Ao = 7 + Cdkl (9) 

with Cd = limA„^o EqX(q)V(l 7 Aox(q))'. This 
constant is generally a decreasing function of the spatial 
dimension d and approaches zero in the limit d —> oo via 
Cd ^ W'^ /8d. Due to the existence of the diffusion pole 
in Eq. (|6b|) the constant Cd becomes infinite in d < 4. 
Equation ^ derived from a Taylor expansion in the local 
vertex Aq does not survive in this form to low dimensions. 

Equation © has generally three solutions for Aq {E + 
iO^,E — iO+). For sufficiently small disorder strengths, 
7 < 7c, all three solutions are real. A perturbative solu- 
tion is of order 7, while two nonperturbative solutions are 
of order ±^/T/Cd- The perturbative solution increases 
and the module of the nonperturbative ones decreases 
with increasing the disorder strength. At a critical ran- 
domness SCrfAg = 1, or equivalently 7c = ^4/27Cd, the 



two positive solutions merge and move into the complex 
plane for 7 > 7c. Disappearance of positive solutions for 
Ao(£' -I- iO+, E — iO+) leads to suppression of the diffu- 
sion pole and simultaneously to vanishing of the diffusion 
constant. Quantity QAq{E + iO~^, E— iO+), emerging be- 
yond the critical point in the localized phase (7 > 7c), 
plays the role of an order parameter for Anderson local- 
ization, see Fig.^ A typical phase diagram for localizcd- 
extended states calculated from Eq. ^ is plotted in 
Fig. 121 Although the mean-field equation © does not 
predict the precise position of the mobility edges, it de- 
termines the mean-field universal properties accurately. 

To conclude, we derived a mean-field approximation 
for two-particle irreducible vertices motivated and justi- 
fied by the asymptotic limit to high dimensions. We suc- 
ceeded in deriving an algebraic equation for the local ir- 
reducible vertex with a bifurcation point at which the dif- 
fusion constant vanishes and a real irreducible vertex be- 
comes complex. A fully consistent and controllable mean- 
field-like theory of the disorder-driven vanishing of diffu- 
sion and Anderson localization was thereby achieved. It 
correctly reproduces the low and high-dimensional limits 
and allows for further systematic improvements. 
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